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Abstract 

We study fluctuations in the drag force experienced by an object moving 
through a granular medium. The successive formation and collapse of jammed 
states give a stick-slip nature to the fluctuations which are periodic at small 
depths but become "stepped" at large depths, a transition which we interpret 
as a consequence of the long-range nature of the force chains and the finite 
size of our experiment. Another important finding is that the mean force and 
the fluctuations appear to be independent of the properties of the contact 
surface between the grains and the dragged object. These results imply that 
the drag force originates in the bulk properties of the granular sample. 
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I. INTRODUCTION. 



The propagation of stress in granular media presents a complex problem which has been 
the subject of extensive study. An applied external stress results in the development of 
an internal structure resisting the stress, a so-called "jammed state" that is dependent on 
the direction and magnitude of the stress The origins of this jamming lie in the fact 

that the forces do not propagate uniformly through the granular sample but are localized 
along directional force chains P ^IO|] , and the jammed state is ultimately characterized by 
the properties of the network of these chains. 

In this paper we investigate the jamming of spherical granular media by analyzing the 
drag force opposing the movement of a solid object through such media. The object's motion 
is opposed by jamming of the grains in front of the object, but if the applied force exceeds 
a critical threshold corresponding to the strength of the jammed state, the object moves 
ahead while displacing the surrounding grains. Such successive breakdowns result in large 
fluctuations in the drag force at low velocities which reflect the strength of successive jammed 
states. Our results demonstrate that the forces arising from such jamming of the grains and 
subsequent bulk grain reorganization dominate the drag process while the frictional forces 
at the surface have little effect. Furthermore the nature of the fluctuations, and thus the 
strength of the jammed states, can be affected by the finite size of the containing vessel 
due to the long-range nature of the force chains. Some of the results have been published 
previously |TT] . 



II. DESCRIPTION OF THE APPARATUS 

The experimental apparatus, shown in Fig. |l], consists of a vertical steel cylinder of 
diameter dc inserted to a depth if in a cylindrical container filled with glass spheres. A 
comb made of three 5 mm diameter steel rods separated by 20 mm gaps is also inserted at 
150 mm depth opposite the vertical cylinder with the role of randomizing the medium. 

The container rotates with constant angular speed while the vertical cylinder is attached 
to an arm that may rotate freely around the rotational axis of the container. Opposing the 
cylinder we have a fixed precision force cell [0 which measures the force F{t) acting on the 
cylinder as a function of time. We also incorporate a spring of known spring constant, fc, 
between the cylinder and the force cell with k varying between 5 to 100 N/cm. The purpose 
of the spring is twofold, it allows the force to slowly build up, and, with a suitable choice of 
k it dominates the elastic response of the cylinder and all other parts of the apparatus so 
that the nonlinear deformations will not significantly affect the results. We vary the speed 
[v) from 0.04 to 1.4 mm/s, the depth of insertion [H) from 20 to 190 mm, and the cylinder 
diameter (dc) from 8 to 24 mm, studying glass (p = 2.5 g/cnr") spheres of diameter {dg) 0.3, 
0.5, 0.7, 0.9, 1.1, 1.4, 1.6, 3.0 and 5.0mm. The force is sampled at 150 Hz and the response 
times of the force cell and the amplifier are less than 0.2 ms. 

As we will discuss in detail later, the drag force experienced by the vertical cylinder, 
F(it), is not constant, but has large stick-slip fluctuations consisting of linear rises associated 
with a compression of the spring and sharp drops (sudden decompressions of the spring) 
corresponding to the collapse of the jammed grains opposing the motion. In most of our 
measurements the angular deflection of the arm during the fluctuations was around 1° (~ 
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one bead diameter), therefore in subsequent discussions we have approximated the spatial 
deflection of the cylinder as linear. Note that these experiments are conducted in dense static 
granular media as opposed to drag in dilute or fluidized media which have been studied both 
theoretically JTSf and experimentally fl^. In our experiments, however, all grains are at rest 
relative to each other and the container for the vast majority of the time. 

During the stick process the cylinder is stationary relative to the granular medium moving 
with speed v. During this process the spring is compressed at a uniform rate. Thus the 
change in the force on the cylinder is described by AF = kAx = kviguck- The slip processes, 
during which the cylinder moves relative to the medium, are very short when compared to 
the duration of the stick processes tghp <C tguck for all depths. Although tguck changes with 
depth {AF has H dependence) tsHp stays approximately the same, and it is comparable 



to the characteristic inertial time = 2t:^ M/k ^ 0.1s. As long as tsuck is significantly 
larger than tgUp the stick slip process probes the properties of the material but if k is chosen 
very large so that the quantity tsUck = AF/kv is comparable to Tin the dynamics will be 
dominated by the inertial properties of the apparatus. Additionally, k must be chosen such 
that it is much smaller than the elasticity of the apparatus. These conditions impose certain 
limits for acceptable values of v and k in a somewhat similar way as with sliding surface 
friction ||T^. As we will show later, however, within these limits the average force and 



fluctuations are not affected by the choice of either v or k. 



III. DATA ANALYSIS 

During each fluctuation the force first rises to a local maximum (the stick process) then 
drops sharply as the jammed state collapses (the slip process). As shown in Fig. |^ we have 
identified three characteristic types of stick-slip motion. These three types, as explained in 
detail below, are: 

• Periodic: The signal resembles an ideal sawtooth signal. The distributions of maxima 
and minima have narrow widths that do not overlap within two standard deviations. 
We observe this behavior for beads of diameter dg = 0.3, 0.5, 0.7, 0.9 and 1.1 mm and 
depths less than 100mm. 

• Random: The signal resembles a random sawtooth signal. The distributions of max- 
ima and minima overlap within one standard deviation. Only beads of diameter larger 
than 1.1 mm, dg = 1.4, 1.6, 3.0 and 5.0 mm exhibit this behavior. We attribute the 
appearance of random fluctuations to the finite size effects associated with both the 
container size and the cylinder/bead ratio. Because of the complications introduced 
by these effects, in the present paper we focus on the other two regimes and we plan 
to study this random regime in more detail in another experiment. 

• Stepped: The force builds up as a sum of small sawtooth like steps. Each subsequent 
maximum and minimum tends to be at a higher force than the previous one, until a 
large reorganization occurs that resets the system and the incremental buildup starts 
again. This behavior can be observed at high depth, {H > 140 mm), and the transition 
from the periodic to the stepped regime can be well characterized. As we will show 
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later, the container diameter has an influence on the critical depth at which the periodic 
regime changes into the stepped one. 



In each of these regimes the quantity in the stick processes stays constant and equal 
to k, thus serving as strong evidence that during the stick processes the cylinder is at rest 
relative to the granular reference frame. We note that we did not observe precursors to the 
slip events even for the lowest speeds (0.05 mm/s) and weakest springs (5 N/m), cases that 
offered the best time resolution. We believe that the gradual, plastic yield observed right 
before slipping of 2D sheared granular layers |]18| does not take place in our system. We 
have extensively studied both the periodic regime and the transition from the periodic to 
the stepped regime. 



A. The Periodic Regime 

1. General characteristics. This regime has been observed for all beads with dg < 1.4 
mm and depths H smaller than 80 mm. A characteristic dataset and power spectrum (the 
squared amplitudes of the Fourier components of F{t)) is shown in Fig. ^ 

As seen on the spectrum, the low frequency interval is dominated by the Fourier terms 
corresponding to the periodicity of the signal. We observe only one or two higher order 
harmonics, while at higher frequencies the spectra exhibit power law scaling with exponent 
of -2. Similar scaling was also observed in sheared granular materials in a Couette geometry 
^ , and is intrinsic to random sawtooth signals . At high frequencies we see only electronic 
noise. 

In the periodic regime, the signal F{t) is made up by continuous rises to a maximum 
value (Fmax) then drops to a minimum (-F^m)- As shown on Fig. ^, the widths of the 
distributions for the maxima and minima are roughly the same and the histograms do not 
intersect. We may consider the separation between the mean value of maxima and minima 
as a parameter in distinguishing the periodic regime from the random one. This parameter 
has a purely statistical interpretation, 4cr separation meaning that ~96% of the minima were 
smaller than f«96% of the maxima. The size distribution of upward and downward jumps 
(the difference between the value of consecutive minima to maxima and maxima to minima) 
is quite similar, as shown in Fig. ^. 

We observed no significant correlation between consecutive values of maxima and min- 
ima, thus we may consider them independent random variables. The interval-averaged val- 
ues, computed as the average value of the maxima or minima over an interval (50 points), 
however, show a correlated overall modulation with an amplitude of about ±5% of the total 
value. Since this modulation affects identically both the maxima and the minima, it appears 
to be caused by small changes of the jamming properties of the media in different parts of 
the container. Such changes might be caused by inhomogeneities in the packing or surface 
level associated with inserting and removing the cylinder. 

2. Depth, cylinder diameter and grain diameter dependence. In a previous paper |17] we 
have reported that the drag force on a cylinder with diameter dc inserted to a depth H in 
a granular bed can be written as F = rjpgdcH'^ where r/ characterizes the grain properties 
(surface friction, packing fraction, etc.), p is the density of the glass beads, and g is gravita- 
tional acceleration. We verified this expected dependence of F on dc and H as shown Fig. 
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p. The power spectra as function of both depth and cyhnder diameter are also presented in 
Fig. 0. 

When we varied the grain diameter, dg, our results were always in agreement with our 
general formula: F = ripgd^H'^. In this formula rj only describes the surface properties 
of the granular material, therefore, the average force is expected to be independent of the 
grain size dg. As we show in Fig. ^. while there is some scattering around the mean for the 
experimentally measured values of rj there is no trend corresponding to the change in grain 
diameter. The period, the time interval between two peaks, is connected to the change in 
the force via the formula T = AF/kv and also shows a linear dependence on cylinder size 
(Fig. 1^) and a quadratic dependence on depth up to /i = 80 mm. 

3. Scaling with k and v. Since our apparatus allowed us to vary k and v, we first verified 
whether the overall behavior of the medium is influenced by the choice of these parameters. 
From a methodological standpoint, decreasing speed or decreasing k affect the results in a 
similar way, by reducing the number of stick-slip events within a certain time interval, thus 
increasing the time resolution of all such events. We were able to vary k between 5 N/cm 
to 800 N/cm, and we find that above k = 100 N/cm the signal started to lose its shape and 
became more and more distorted with increasing k (see Fig. plQl) . 

Since the amplitude of the signal and thus the period is inversely proportional to k, 
increasing k reduced the duration of the process during which the cylinder moved with the 
media {tguck) to the point at which the inertial effects and the elasticity of the jammed 
medium and apparatus were interfering with the force measurements. This has set an upper 
limit for acceptable values of k. Below this limit, however, we could not observe any influence 
of k on the data. We verify this assumption by rescaling the spectra corresponding to the 
formula AF = kvAt. If k and v affects the data only via this relation then the scaled power 
spectra expressed as kvPiu) = f{u/kv) should overlap. 

As seen in Fig. |lT], the power spectra collapse very well with k (for data taken with 
constant v). A similar scaling performed for the speed v, with v varied from 0.05 mm/s 
to 0.5 mm/s is shown in Fig. |l^ {k fixed). Based on these two results we believe that 
the choices of k and v did not affect the nature of the force fluctuations measured by this 
experiment, i.e. that the force fluctuations represent static properties of the grains rather 
than the measurement process. 

Surface friction and cylinder profile. We have varied the friction coefficient of the surface 
of the cylinders by using cylinders made from different materials. We tested sandblasted 
aluminum and steel, smooth plexiglass and teflon coated steel cylinders. By this substitution 
we varied the static friction coefficient by a factor of ~ 2.5 The mean drag force was 
almost identical for all of these cylinders (within 5%) and the fluctuations also showed no 
qualitative difference. (Fig. p!3|). 

Thus, it appears that the frictional forces at the surface of an object have a negligible 
contribution to the drag force. This implies that the fluctuations originate in the bulk of 
the medium. It is also noteworthy that the mass of these cylinders varied significantly (the 
steel was around 3 times heavier than the aluminum cylinder) demonstrating that inertial 
effects are not significant. 

We have also tested the effects of the shape of the object on the drag force, and we 
found similarly unexpected behavior. There was very little difference in the signal profile 
when we substituted our circular cylinder with a half cylinder with the plane of the cut 
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perpendicular to the direction of motion. While we do not want to suggest that the force is 
totally independent of the shape of the object we believe that the shape has a much smaller 
impact than one would expect from an analogy with the shape effects on the drag in fluids 
(Fig. HI). 

This relative independence of cylinder shape and surface properties can, however, be 
explained quahtatively as a consequence of the nature of force propagation in granular 
media. As the cylinder advances it has to displace material through dilation, thus it has to 
reorganize the grains into spatially different positions. Since the force chains have a long- 
range character, the volume of the material resisting the movement is much larger than the 
volume of material engaged into frictional interactions at the surface of the object. The 
intergrain interactions are not only independent of the cylinder-grain interactions but they 
are also far more numerous. Thus, it is easy to see that they may completely dominate the 
dynamics. One can also describe this mechanism based on the energetics of the stick-slip 
process. The energy stored in the spring will be dissipated in friction at the surface of the 
object and in friction among grains far from the cylinder, as the work necessary to bring 
them into their new spatial positions. The number of the grains involved is not affected 
by the friction coefficient of the object and is only slightly influenced by the shape of the 
object. 



B. Stepped regime 

A striking feature of the data is that the fluctuations change character with depth. For 
smaller grains dg < 1.6 mm and H < ^ 80 mm the fluctuations are quite periodic, i.e. 
F{t) increases continuously to a nearly constant value of F^ax and then collapses with a 
nearly constant drop of AF (Fig. |T5|) as discussed above. As the depth increases, however, 
we observe a change in F{t) to a "stepped" signal: instead of a long linear increase followed 
by a roughly equal sudden drop, F{t) rises in small linear increments to increasing values of 
Fmaxi followed by small drops (in which AF is significantly smaller than the rises, 15-20%), 
until Fjnax reaches a characteristic high value, at which point a large drop is observed. This 
transition from a periodic to a "stepped" regime is best quantified in Fig. |16|, where we plot 
the depth dependence of the mean of AF. In the periodic regime, AF increases with depth 
just as F does. As the large uniform rises of the periodic regime are broken up by the small 
intermediate drops, however, AF shows a local minimum and then increases continuously 
again. The transition can also be observed in the power spectra as shown in Fig. |^ For low 
depths the power spectra display a distinct peak characteristic of periodic fluctuations, but 
these peaks are suppressed for H > He in correlation with the changes in the qualitative 
character of F{t). 

The transition can also be observed in the histograms (curves that can be well approxi- 
mated as gaussians) depicting the distribution of maxima and minima (Fig. |1^). The initially 
well-separated peaks corresponding to an ideal sawtooth pattern widen while the distance 
between their maxima decreases up to the point at which they intersect within less than one 
a separation. 

Another useful way of characterizing the structure of the fluctuations and the transition 
to the stepped regime is by building return maps, plotting F^"^ vs F^"f (Fig. |l^). In these 
plots, the continuous line shows the F^"^ = F^^f line, points above it mean that the A^+lth 
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maximum is larger than the A^th one, while points below mean that the A^th maximum is 
larger than the next one. 

As long as there is no relation between consecutive maxima and minima these maps 
should not show any characteristic shape and should form a circular pattern with gaussian 
spread. Indeed at low depths we do observe such a pattern but with a somewhat elongated 
structure, the elongation being probably caused by the previously mentioned overall modu- 
lation of the signal. At high depths, however, the blob starts to show a structure associated 
with the stepped character of the motion, suggesting a history dependence in the way each 
jammed state will further evolve. 

The graphs show that at large depths the majority of the peaks are followed by a peak 
with an even greater height, and the increase from one peak to the next one is within an 
interval that has a well defined upper limit. Interestingly this limit forms a line shifted above 
and parallel with the line of F^""^ = FJ^!^f. Only at the highest values of -Fmax do we observe 
values that are smaller, usually significantly smaller than their predecessor. This means that 
the force always has to first build up to a certain value before a large reorganization can 
occur. 



C. The transition from periodic to stepped fluctuations. 

The transition from a periodic to a " stepped" signal is rather unexpected, since it implies 
qualitative changes in the failure and reorganization process as H increases and the existence 
of a critical depth. He- An explanation for He could be provided by Janssen's law ||20| which 
states that the average pressure (which correlates directly with the local failure process) 
should become depth independent below some critical depth in containers with finite width. 
This should not occur in our container, however, which has a diameter of 25 cm, significantly 
larger than He (— 10cm). Furthermore we see no deviation in the behavior of F{H) from 
F oc H^, which depends on the pressure increasing linearly with the depth (Fig. |l^ inset). 

In order to account for the observed transition, we must inspect how the force chains 
originating at the surface of the cylinder nucleate the reorganizations. The motion of the 
cylinder attempting to advance relative to the grains is opposed by force chains that start 
at the cylinder's surface and propagate on average in the direction of the cylinder's motion. 
These force chains will terminate rather differently depending on the depth at which they 



originate, as shown schematically in Fig. |19[ For small H, some force chains will terminate 
on the top surface of the granular sample and the stress can be relieved by a slight distortion 
of the surface. Force chains originating at large depths, however, will all terminate at the 
container's walls. Since the wall does not permit stress relaxation, the grains in these force 
chains will be more rigidly held in place. According to this picture. He corresponds to 
the smallest depth for which all force chains terminate on the wall. When the cylinder 
applies stress on the medium, the force chains originating at small H {H < He) reduce their 
strain through a microscopic upward relaxation of chains ending on the free surface. By 
contrast, the higher rigidity of force chains originating aX H > He impedes such microscopic 
relaxations. Thus a higher proportion of the total force applied by the cylinder will be 
supported by those force chains, enhancing the probability of a local slip event occurring at 
high depths. Such a slip event would not necessarily reorganize the grains at all depths (for 
example the grains closer to the surface may not be near the threshold of reorganization). 
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thus the shp event might induce only a local reorganization and a small drop in F{t). The 
large drops in F{t) would occur when force chains at all depths are strained to the point 
where the local forces are close to the threshold for a slip event. This scenario also explains 
why F{H) does not change at He, since F is determined by the collective collapse of the 
jammed structure of the system. 

According to this picture, the transition is expected at smaller depths in smaller con- 
tainers since the force chains would terminate on the walls sooner (seen Fig. |T^). Indeed, 
as we show in Fig. 0, the transition does occur at a depth approximately 20 mm smaller 
when the measurements are performed in a container 2.5 times smaller (with diameter of 
100 mm). Furthermore, we fail to observe the periodic fluctuations in any grains with di- 
ameters 1.4 mm or larger, which is consistent with the suggested mechanism, since larger 
grains correspond to a smaller effective system size. 

In Fig. ^ we show how for low depths the signal changes from nearly periodic to non- 
periodic just by changing the grain diameter. Similarly, the high depth regimes that show 



a stepped behavior become periodic for grains of smaller diameter (Fig. ^T]). This result 



is consistent with our earlier observation that reducing the container size has resulted in a 
decrease of the He limit. 



IV. CONCLUSIONS 

In summary, we have made a detailed study of the fluctuations in granular drag. We find 
that the fluctuations are strongly affected by the long-range nature of the force propagation 
in granular media. The jamming responsible for these fluctuations originates from a localized 
applied stress and our results point towards the need of better understanding of how force 
chains originating from a point source propagate and disperse geometrically. 

We gratefully acknowledge the support of the Petroleum Research Foundation adminis- 
tered by the ACS, the Alfred P. Sloan Foundation, NSF grants PHYS95-31383 and DMR97- 
01998 and NASA grant NAG3-2384. 



8 



REFERENCES 



e-mail: schiffer@phys.psu.edu; 
[1] M. E. Gates, J. P. Wittmer, J.-P. Bouchaud and P. Claudin, Phys. Rev. Lett. 81, 1841 

(1998) and |cond-mat/9901009| (1999); M. E. Gates and J. P. Wittmer, Physica A 263, 
354 (1999). 

A. J. Liu and S. R. Nagel, Nature (London) 396, 21 (1998). 

C. Liu S. R. Nagel, D. A. Schecter, S. N. Coppersmith, S. Majumdar, O. Narayan and 
T. A. Witten, Science 269, 513 (1995); D.M. Mueth, H. M. Jaeger, and S. R. Nagel, 
Phy. Rev. E. 57, 3164 (1998). 

B. Miller, C. O'Hern, and R. P. Behringer, Phys. Rev. Lett 77, 3110 (1996). 
X.Jia, C. Caroli, B. Velicky, Phys Rev. Lett. 82, 1863 (1999). 
E. Kolb, T. Mazozi, E. Clement, and J. Duran Eur. Phys. J. B. 8, 483 (1999). 
S. N. Coppersmith C.-h. Liu, S. Majumdar, O. Narayan and T. A. Witten, Phys. Rev. 
E 53, 4673 (1996); M. L. Nguyen and S. N. Coppersmith 59, 5870 (1999). 
A. V. Tkachenko and T. Witten, |cond-mat/991025q . 

L. Vanel, D. Howell, D. Clark, R. P. Behringer and E. Clment, Phys. Rev. E. 60, 5040 

(1999) ; D. Howell, R. P. Behringer and C. Veje, Phys. Rev. Lett. 82, 5241 (1999). 
K. Wieghardt, Annu. Rev. Fluid. Mech. 7, 89 (1975). 

I. Albert, P. Tegzes, B. Kahng, R. Albert, J. G. Sample, M. Pfeifer, A.-L. Barabasi, T. 
Vicsek and P. Schiffer, Phys. Rev. Lett. 84, 5122 (2000) 
Sensotec, Precision Miniature Load Cell, Model 31, 10 lbs range. 
V. Buchholtz, T. Poschel, Gran. Mat. 1, 33 (1998). 
O. Zik, J. Stavans, and Y. Rabin Europhy. Lett. 17, 315 (1992). 

It is an experimental observation of sliding friction that stick-slip behavior disappears 
if the spring k is stiff enough, or if the sliding velocity v is high enough. For an overview 
article see B. N. J. Persson, Surf. Sci. Rep. 33, 33 (1999). 

A. L. Demirel and S. Granick, Phys. Rev. Lett. 77, 4330 (1996); H. J. S. Feder and J. 
Feder, Phys. Rev. Lett. 66, 2669 (1991). 

R. Albert, M. A. Pfeifer, A.-L. Barabasi, and P. Schiffer, Phys. Rev. Lett. 82, 205 
(1999). 

S. Nasuno, A. Kudrolh, and J. P. Gollub, Phys. Rev. Lett. 79, 949 (1997); S. Nasuno, 
A. Kudrolli, A. Bak, and J. P. Gollub, Phys. Rev. E 58, 2161 (1998); J.-C. Geminard, 
W. Losert, and J. P. Gollub, Phys Rev. E 59, 5881 (1999). 

We measured the friction angle (p for the cylinders sliding on a plate covered with glass 
beads glued onto its surface. (/i = taii{ip)). 

R. L. Brown and J. C. Richards, Principles of Powder Mechanics (Pergamon Press, 
Oxford, 1970). 

This is why we chose to explain the periodic - stepped transition in the context of 
container size and not grain diameter. 

A. Ngadi amd J. Rajchenbach, Phys. Rev. Lett. 80, 273 (1998) 



9 



FIGURES 



FIG. 1. Schematic view of the apparatus used for measuring the drag force. The grains were 
contained within a 25 cm diameter rotating bucket as described in detail previously |]ll| a) top 
view b) cross-sectional view. 

FIG. 2. Graphical representation of the drag force fluctuations in the three observed regimes, 
periodic, random and stepped. Transitions between these regimes are determined by parameters 
like grain size, insertion depth and cylinder diameter. 

FIG. 3. The signal shape and the power spectrum in the periodic regime {dc = 10 mm , dg = 0.9 
mm and depth H = 80 mm) . The lower graph shows a typical power spectrum for this regime. 
The dashed line has the slope of —2. 

FIG. 4. Distribution of maxima and minima in the periodic regime, (dc = 16 mm , dg = 0.9 
mm and H = 60 mm) . The large separation between the peaks of the histograms for the maxima 
and minima characterizes the periodic regime. 

FIG. 5. Fluctuation (upward rise and downward step) distribution in the periodic regimes for 
dc = 16 mm , dg = 0.9 mm and depth H = 80 mm . The upward rises correspond to height 
differences from a minimum to the next maximum. Downward steps correspond to changes from 
a maximum to the next minimum (absolute values). 

FIG. 6. The mean value of the drag force as a function of diameter and depth. The force has 
a quadratic dependence on the depth and a linear dependence on the diameter of the cylinder 
{dg = 0.9 mm) . 

FIG. 7. Power spectra for different depths with cylinder of diameter dc=16 mm (upper graph) 
and different diameters at depth H = 80 mm (lower graph), dg=0.9mm. Data are vertically offset 
for clarity 

FIG. 8. The drag coefficient, r/ =< F/dc > /pgH"^, as function of bead diameter. These data 
are for H = 80, and for each datapoint we averaged over 5 cylinder diameters. The error bars 
represent the standard deviation. 

FIG. 9. Dependence of the size of upward rises (minima to next maxima) in the force on the 
cylinder diameter (dg=0.9 mm, H =80 mm). 

FIG. 10. The effect of increasing k on the signal. At very high spring stiffness the nonuniform 
elasticity of the apparatus will dominate the response and the inertial effects play an increasing 
role H = 170 mm , dc = 10 mm , dg = 0.9 mm . 
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FIG. 11. Scaling with k in the periodic regime. The spring constant was varied from 25N/cm 
to lOON/cm with the leftmost peak corresponding to the weakest spring {H = 80 mm , c?c = 16 
mm , dg = 0.9 mm) . 

FIG. 12. Scaling with v in the periodic regime. The speed was varied from 0.05 to 0.2 mm/s 
with the leftmost peak corresponding to the lowest speed {H = 60 mm , dc = 16 mm , dg = 0.9 
mm) . 

FIG. 13. Effects of the surface friction of the cylinder on the mean force and fluctuations 
{dc = 16 mm , dg = 0.9 mm). The error bars represent the standard deviations, 

FIG. 14. Comparison of the average force for different surface friction and shape {dg = 0.9 mm) . 

FIG. 15. The characteristic fluctuations in the drag force at 4 different values of H for dc = 10 
mm. Note the transition from purely periodic fluctuations {H < 60) to stepped fluctuations with 
increasing depth (if > 100). 

FIG. 16. The transition from periodic to stepped fluctuations as shown through the magnitude 
of the average drop AF, for two different container sizes: circles - large container {D =25mm), 
triangles - small container {D = 10 mm). The transition occurs at smaller He in the smaller 
container {dc = 10 mm). The inset shows the depth dependence of the average drag force for 
dg = 1.1 mm and dc = 10 mm, and there is no change in slope at He- The solid line has slope 2. 

FIG. 17. Histograms for low depth and high depth regimes. The periodic regime shows well 
separated peaks while the distributions in the stepped regime overlap within one standard devia- 
tion. 

FIG. 18. The i^jv"i vs Fjy"^ return maps for dg = 0.9 mm , dc=10 mm and different depths. 
The straight line is the F^"^ = FJ^'^^ line. Points above the line indicate an increasing trend of 
the rises, the points below show a decrease. 

FIG. 19. Schematic diagram of the force chains perturbed by the presence of the walls. While 
the figure is greatly simplified since the force chains bifurcate and follow nonlinear paths there is 
every reason to believe that their effective area has a conical shape. 

FIG. 20. Low depth {H < He = 100 mm ) signal shape for different grain diameters. The small 
grain size show a well defined periodicity while the large grain sizes manifest a random pattern 
{dc = 16 mm) . 
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FIG. 21. High depth (H > He = 100 mm ) signal shape for different grain diameters. The 
changes due to the increase in grain size cause the pattern to turn from nearly periodic to the 
stepped stick-slip motion {dc = 10 mm , dg = 0.9 mm) . 
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